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We study matrix models related via the correspondence of Dijkgraaf and Vafa to supersym¬ 
metric gauge theories with matter in the fundamental. As in flavorless examples, measure 
factors of the matrix integral reproduce information about R-symmetry violation in the 
held theory. The models, studied previously as models of open strings, exhibit a large-M 
phase transition as the number of havors is varied. This is the matrix model’s manifes¬ 
tation of the end of asymptotic freedom. Using the relation to a quiver gauge theory, 
we extract the effective glueball superpotential and Seiberg-Witten curve from the matrix 
model. 
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1. Introduction and summary 

Dijkgraaf and Vafa have found that matrix integrals compute all of the holomorphic 
information in A/” = 1 gauge theories with classical gauge groups 0,i|,|l] • In the interest of 
understanding their proposal better, we will look at the matrix model that they prescribe 
for M = 2 gauge theory with Nf fundamental hypers mass-deformed to A/" = 1 by a tree 
level superpotential for the adjoint chiral held $, and masses for the squarks. These held 
theories, studied in Q, exhibit fascinating behavior as a function of the number of havors. 

Such models arise in string theory in a number of ways. The U{N) theories arise by 
probing the “canonical example” of Dijkgraaf-Vafa with D5-branes on a noncompact curve 
of the Calabi-Yau (CY) [§. Flavorful theories with real gauge groups arise on N D3-brane 
probes of D7-brane conhgurations |P,|^,||,^; this fact was recently employed to add holes 
to the BMN worldsheet H- We will instead use it to poke holes in the random surfaces 
described by the Dijkgraaf-Vafa matrix integrals. 

In particular, for N = 1, the theories with Nf < 4 (with Sp{2N) = SU{2) gauge 
group) arise on a D3-brane probe of a resolved singularity of F-theory. The theory 
with Nf = 4 massless flavors is obtained by creating the singularity; in this case the 
dilaton is constant and the gauge theory is scale-invariant. This brings us to the question 
we would like to answer: How is this dependence on Nf manifested in the Dijkgraaf-Vafa 
matrix model? 

We are therefore led to consider a matrix integral of the form 

I (1.1) 


Z{gk,m^,M) = / dMQdQexp -Wo{4>) + 


Mf 

a=l 


where 

n-\-l 

Wom = J29ktr 

fc=i 

$ is a complex M x M matrix, Q is a complex M x Mf rectangle, and Q is a complex 
Mf X M rectangle. An arbitrary meson source can be brought into this diagonal 

form by rotating the Qs. As in the work of Dijkgraaf and Vafa, these are to be thought of 
as line integrals over matrices of complex numbers. 

Without using very much technology, we can study in detail the model which arises 
upon integrating out the fundamentals. This generates a ln(m — $) potential for $. The 


model with this addition to the potential can still be solved by the method of BIPZ |11 
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and actually was studied in this way as a model for open strings by Kazakov [1^|0. The 
coefficient of the log is 7 = and we will show that the cuts close up when 7 is chosen 

so that the gauge theory can be conformal. Related by supersymmetry to this statement 
is the fact that the R-symmetry becomes non-anomalous for this choice of Nf/N. This 
fact can be detected in the matrix model through the dependence of the measure on 7 . 
Finally, we extract the glueball superpotential and Seiberg-Witten curve from the large-M 
solution to the matrix model. 


M’s and N’s 


One of the more mysterious aspects of the Dijkgraaf-Vafa prescription, at least to the 
author, is the disjunction between the number of colors N of the gauge theory, and the 
number of colors M of the matrix model (which plays the role of the glueball superfield). 
The addition oi Nf flavors to the gauge theory only complicates this issue. We will in¬ 
troduce a number of flavors M/ in the matrix model which is again not the same as the 
corresponding number in the gauge theory. We will, however, identify the ratio 

^ N M ’ 

This will be the parameter of interest in our study of the matrix model. We take this as part 
of the prescription, but (thinking of 7 as the weight with which holes in the random surface 
contribute) one which is again motivated by topological string duality 0 . Specifically, 
this is the usual Dijkgraaf-Vafa limit for a two-node quiver gauge theory [H which reduces 
to the flavorful theory when the dynamics of the gauged flavor group are frozen out. It is 
from this perspective that we will be able to extract the superpotential. 


Related work 


Matrix models with fundamentals in this context are mentioned in a footnote in 


They also make an appearance in the very recent |17,18|. Other work on understanding 


and extending the Dijkgraaf-Vafa proposal includes ||19|,^|^,|2^,^,|24|,| 


^ Related models were also studied in |13], [14|, and in particular in [15| where the critical 
behavior at 7 = 2 was explored. 
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2. Flavorful matrix models 


To begin, we write down the matrix integral with naive conplings and helds marked 
with hats: 


Mf 


Z= dddQdQexp -lTo($) + 


rrir 


( 2 . 1 ) 


a=l 


with 


ITo(^) = .^itr $ + ^2tr + .... 


q; is a flavor index; color indices will be denoted a,b,... These hatted flelds and conplings 
will be related to those which shonld have flnite large-M limits (which lack hats) by an 
M-dependent rescaling. These hatless variables are chosen below so that there is a well- 
peaked saddle point of the $ integral, the location of which is M-independent. We will 
observe that the precions 

-M^ InM 
2 

term in the matrix model free energy, which is derived from the inverse volnme of the 
matrix model gange gronp, can also be detected by snch a propitions held rescaling. 


2.1. Integrating out the flavor 


For the moment, we are interested in the regime of couplings where the quark masses, 
nia, are much bigger than the bare mass of the adjoint in kFo(fl’). In this regime, we 
integrate out the fast Q modes at fixed $ to get an effective potential for fl>. The integrals 
over Qa are Mf independent gaussian integrals. This gives 




a.=l 

( ^ ^ ^ \ 
(i$exp j —lTo($) — ^ trln($ — rhal) 

\ / 


Setting all of the masses equal to m for simplicity, this is 




(i$ exp ( —IFo(fl>) — M/tr ln($ — m) 


( 2 . 2 ) 


(2.3) 


A matrix integral very similar to (|2.3| ) was studied by |T1,|1,|T1,0 as a discretization 
of an open string worldsheet. In this model, the counterpart of M/ is the weight accom¬ 
panying a hole insertion. The logarithmic potential was chosen to reproduce a sum over 
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discretizations of the worldsheet boundaries, with equal weight for arbitrary numbers of 
segments of the boundary. 

Other than a relabeling of couplings, the difference between our model and that of 
Kazakov is that the logarithmic potential term of is 

tr ln(m — (p^). 

This is the potential that would arise if the Af = 2 superpotential were instead of 

Q^Q. The held redehnition $ = 99 ^ required to relate the two integrals directly introduces 
a Jacobian factor which adds a term 


1 


-tr ln$ 


to the potential. From the calculation ( p.2|) above we see that this is the same as the 
effect of adding M/2 massless hypers. We will hnd it convenient to solve the integral (|2.3| ) 
directly. The qualitative behavior we hnd is the same as that found in . 

2.2. The continuum 

In order to proceed, diagonalize the matrix $ as $ = WDW with 

T> = diag(Ai,...,AM), 

and U unitary. These eigenvalues are normalized as in [|n|. By the magic of logarithms, the 


integrand of (|27^ ) does not depend on the angular U variables. Their integration produces 
the Vandermonde determinant 


AiX) = l[{Xa-Xtf 

a<b 


The integral becomes 


. Ml Mf 

Z = / JJdA„A(A)exp^ I -lTo(A„) + ^ ln(A„ - 

a a=l \ a=l 


(2.4) 


At this point, it is convenient to introduce a continuum in the space of colors. Let 

1 ^ 

Xa = VMX{a = a/M), da = — . 

a=l 


The eigenvalue density 


is normalized to 


da 

= lx 


dnp{n) = 1 . 


(2.5) 
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Large M scaling 


We now introduce the promised variables without hats: 

92 . 92. . 9k 


S’ sVm' 

thi = y/Mnii 

1 


<9k = 


9k 




( 2 . 6 ) 


$ = Q = 


uy^gl^^ 


Q, Q = 


;Q- 


Here we have hnally introduced the quantity S = 9 sLL^ which is hxed in the large-M limit. 
Note that these rescalings are closely related to those made on dimensional grounds in 
the stringy realization of the gauge theory [^. Plugging these into ( p.l|) , we hnd that the 
resulting A integral is of the form 


Z = / D\{a) exp ( 4 -^ 0 [A] + C[M) 

K 9 S 


(2.7) 


with C(M) independent of A and 


JFo[A] = da db\n{X{a) — X{b)) — S / da Wo(A(a)) — 5'^7 / daln{X{a) — m) 

Jo Jo Jo Jo 

^ J dxjdzp(\)p(z) lii(A - 2 ) - S j d\p(\]Wo{X) - J dXp{\) ln(A - m) 

( 2 . 8 ) 

Here Hd)(A) = giX + ( 72 A^ + .... The crucial feature of (|2.7|) is that JTo[A] is independent 
of M. This is the normalization used by Dijkgraaf and Vafa; the matrix action is 


(2.9) 


— {Wo{^) + Q^Q - mQQ 
9s 


Now we return to the “constant (/-independent term” llTl[| C{M). This held- and 
coupling-independent term was not relevant for previous applications of matrix integrals. 
It is 

C(M) = In M - ^MfM In M 
= (2 -7)^M2lnM. 

This reproduces the leading M-dependence of the log of the inverse volume of U{M) 
in the held normalization we are using. Further, it provides the “entropy factor” arising 
from the havor integral^!. The Dijkgraaf-Vafa prescription relates the matrix model free 


( 2 . 10 ) 


I am grateful to Nissan Itzhaki for comments on this point. 
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energy to the prepotential of the gauge theory. In parallel with the discussion of |]^ this 
term leads to the following contribution to the effective superpotential of the gauge theory: 

ir.//(S) = (2-7)A'Slii^, (2.11) 

up to linear terms in S which are independent of the cutoff scale Aq. The prefactor of 
this Veneziano-Yankielowicz superpotential is proportional to the anomaly in the 
current of the held theory. The introduction oi Nf havors in the fundamental modihes 
this from 2N to (2 — 7 )A^. It is gratifying that this is reproduced by the simple matrix 
integral. 

2.3. Solution at large M 

In terms of the variables normalized to have a hnite large-M limit, the saddle point 
equation is 

Vi(A) + = 2 / (2.12) 

b A — m J A — n 

Rewrite this equation as 

^w;,(\) = 2j (2.13) 

where 

PoifJ^) = pM (2.14) 
or more generally in the case of arbitrary masses 

Poip) = pip) - 7 ^^ 

■f a=l 

Therefore the eigenvalue density po satishes the same integral equation as that of the 
theory without havors, with the modihed boundary condition 

f f f 1 

/ dppoip) = / dpp{p) - dpY^ ~ ^ 

a=l h 

This integral is to be performed over the real p line. 

For simplicity, let us consider the case of a cubic superpotential, in the case of a single 
cut, i.e. choose the potential Wq to have a single critical point. Placing the cut at [2a, 26], 
the BIPZ solution for the resolvent 

u,o(\) = (2.15) 

J2a P ^ 
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of the theory without flavors determines the solution of (p.l2). This is [n 


5'wo(2) = ‘ 2 .g 2 Z + - { 2 g 2 + “igz^a + 6 ) + dfifsz) \/{z - 2a){z - 26). (2.16) 


The conditions determining the positions a, b of the ends of the cut are determined by the 
behavior of (|2.16|) at z —> oo. For the model with flavor, these are 


Sgsib - af + 2(a + b){2g2 + Sgsia + b)) = 0 (2.17) 

^(6 - af{ 2 g 2 + Ggsia + 6 )) = 2 - 7 . (2.18) 

Up to coupling redehnitions, these differ from equations (46) of |]^ only by the replacement 
2 I— > 2 — 7 in the second condition. Note that they do not depend on m. 


Healing of cuts 


Assume there is a stable vacuum at the origin, and place all of the eigenvalues there, 
e.g. consider gk = 0, k > 92 > 0- From (|2.17|) and ( |2.18|) we immediately see that when 

7 —> 2, b — a 0. That is, the cut closes up. Beyond 7 = 2, the cut at the origin moves into 
the complex plane. In the context of a Hermitian matrix integral this was interpreted 
as a large-M phase transition beyond which the theory lacked a stable solution. However, 
when the saddle point we are studying is that of a holomorphic line integral, this is not so 
catastrophic. In fact, the cut surrounding the unstable extrema of W always extend into 
the imaginary direction of A. 

In the matrix models for conhning gauge theories, the size of the cuts which hold the 
eigenvalues goes like the IR scale A of the held theory. The closing of the cuts is a signal 
that the theories become scale invariant, and then no longer asymptotically free as Nf 
passes through 2N. 

We can see the corresponding effect for more general tree-level superpotentials as 
follows. As a consequence of the saddle equation (|2.13|) , the resolvent 

wo(a^) = [ 

J X - X 


of the 7 = 0 theory satishes at large M an algebraic equation of the form \e.g. [28 


uJoixf + -uo{x)Wq{x) -J^foix) = 0. 


S 


4A2' 


(2.19) 
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Here fo{x) is a degree n — 1 polynomial in x, which should be thought of as a function of 
Si = QsMi, with Mi the number of eigenvalues in the zth cut. This polynomial /o differs 
from the one in the 7 = 0 solution through its dependence on the locations of the cuts. 
Note that the resolvent of the theory with flavor is related to ujq by the addition of a pole 
term 

u(x) = ciJo(x) H- - —. ( 2 . 20 ) 

X — m 


The remainder term, /o, can be written [PP 


/o(^) = 4^ 


dw po{w) 


Wi(w)-Wi(z) ^ 

w — z 


this is a polynomial of degree n — 1 in z, fo{z) = The fact that / dw po{w) = 

2 — 7 then implies that as 7 —2, the leading term, bn-i, vanishes. In more detail, (|2.19| ) 
is a quadratic equation for ujq, which therefore has the solution 


2SM^) = w;,{z) ± ywTiTTT). 


From this expression we learn that 


hn-i = -2(2 - 7 )(n + l)gn+iS. 


But, we also know from and from that 


b-n—l 


-4(n + l)gn+i 


dWeff 
d In 


Therefore, we see from the matrix model that the superpotential is independent of the IR 
scale of the gauge theory when 7 —2. If m = 0, all the cuts are healed, in accord with 
the fact that in that case scale invariance is never broken. 


2.4- The Seiberg-Witten curve and glueball superpotential 


In this subsection we will extract the glueball superpotential from the free energy of 
the matrix model, and the Seiberg-Witten (SW) curve of the gauge theory from the loop 
equationl. 


^ Please note that this subsection was rewritten for version three, after the appearance of 
32,32,14,35,36,37,^ in which superpotentials were found using matrix models with a finite 


number of flavors. 































In the canonical example of Dijkgraaf-Vafa, the free energy of the matrix model de¬ 
termines the effective glueball superpotential of the held theory via 


Weff{S,)=Y,N. 


dJ^o 

dS, 


2'KiTS 


( 2 . 21 ) 


with = ^r^lnZ. Thinking of the havors as arising from the M = 2 quiver with 

two nodes of respective ranks N and Nf m. the limit that the coupling of the havor group 
vanishes, this formula is modihed in our case tJ 


W,„ = +Y,N.^ + 


( 2 . 22 ) 


dS^ ^ dSf 

where Sf = QsMf. Because the havor gauge coupling is zero, we do not include the 


dynamics of the corresponding matrix variables, and treat them as a background. (| 2 . 22 | ) 
has been written for the case where all of the havor branes are coincident (he. the masses 
for the havors are identical), but the generalization is clear. 

A convenient way |j^ uw 
to 


to compute is by varying the eigenvalue density according 


p{z) p{z) + SSi^5{z - ei) 


with ei some hxed point inside the zth cut. From 


evaluated in the saddle, we hnd 


dS,. 


fAn 


dx 


Wl){x) - 5 2 dX 


P(A) 
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X — X X — X 


- pW=PoW + ^S{\-m) 


We have ignored the additive constant, and Ao-independent terms linear in S, which can be 
absorbed in the bare coupling, r. This expression is identical to that for the case without 
havor; the contributions proportional to 7 cancel each other. So we have 


dJ^o 

dSi 


= / dxWQ{x) — 2S dwpo{w) 


dx 


'Bi 


ISi X-W 


' Bi 


dx {Wq{x) + 2SijOQ{x)) (2.23) 


where Bi is a contour running from the reference point to the cutoh Aq. As in the case 
without havor, we introduce y, the singular piece of the (7 = 0 ) resolvent, by completing 
the square in ( |2.19|) 0 

y{x) = 2 Sujq{x) + Wq{x). 


^ While this paper was being revised, a closely related observation was made independently in 
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Then (|2.23|) becomes 


9^0 

dS, 


= / ydx. 


'Bi 


Similarly, we can determine the valne of the new term in (p.22|) by varying 

6'y 


p{z) ^ p{z) + y - mo); 


(2.24) 


(2.25) 


mo is the mass at which we add the new flavor. The factor of two in (|2.25|) appears becanse 


of the corresponding factor in (|2.14|) . Under this change, 
1 


= 


-S'lUo(mo) + / dX p{X) ln(A — mo) — ln(mo — m) 


p{\)=po{X) + ^5{\-m} 
(2.26) 


Again up to irrelevant terms, this can be rewritten to give 


dJ^o 1 


dSf 


dxWQ{x) + 2S f dx [ (po{X) + ^S{X-m)]+'yS 
TiQ J rriQ J ^ X \ 2 / 


r'Ao 


dx 


Jmo 

r*Ao 


mo 


X — m 


1 


rAo 


" JrriQ 

We therefore find 


dx{2SijjQ{x) -\-Wq{x)) = I ydx. 


' mo 




r*An 


ydx 


(2.27) 


'Bi 


with y dehned by 


r = W - fo 


Here, the coefficients of /o(ic) = Ylik^k{Si)x^ are determined by the dynamical glueball 
helds Si according to 


= 


ydx 


where Ai is a contour encircling the fth cut. This is the formula for the effective superpo- 
tential predicted from the physics of D5-branes on the generalized conifold. 

Taking = 1, minimization of this superpotential with respect to variations of the 
polynomial /o has been shown |l40| , ^ ^ using the methods of |^l[] to result in the correct 
Seiberg-Witten curve 


N 


y 


= ^^{x-m)^A 


(2.28) 


a=l 


where (pa are the critical points of ITo. 
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3. Discussion and prospects 

In this paper, we have focused on the regime of couplings << mg (though of 
course we can still perform the gaussian integral over Q in the other regime) where the 
interesting observables involve the adjoint held, as in the case without havors. It will 
be interesting to try to compute other observables involving light quarks. In the SU{N) 
version of these theories, these include baryon operators 


Z[M,B,B] = j dMQdQexp f Vko($) - QoQ 


rur 


exp 


[y 13 ^ a + ^ • Vajv + ^ai ...fiiv Vai ' ' ' Vajv j 


These baryon sources exist for Nf > N, corresponding to a value of 7 at which we have 
not yet detected any change in behavior of the matrix model. 


Beyond the transition 

Consider the F-theory realization of the related symplectic models mentioned in the 
introduction {Sp{N) with Nf < 4 and an antisymmetric tensor). The havor symmetry 
(which is SO(8) in the critical case) is the gauge symmetry on the D7-branes. Increasing 
Nf /N beyond the critical value is achieved by adding more D7-branes to the D 4 singularity. 
This is possible without destroying the triviality of the canonical bundle, and one obtains 
in this way collections of D7-branes with the exceptional series of gauge groups of rank up 
to 8 . The D3 probe theories are then held theories with exceptional havor symmetry [e.g. 

@]- 

It has become clear that the complex x-plane of $ eigenvalues can be identihed with 
the image of a hbration of a noncompact CY geometry. A cut which holds the eigenvalues 
in the large M solution is identihed with the image in the x-plane of a three-cycle in 
this geometry (after a geometric transition induced by the hux generating IFo($)). It is 
therefore tempting to speculate that tuning 7 past the critical value is the matrix model 
version of performing an extremal transition in the CY geometry, during which the three- 
cycle shrinks and one hnds an even-dimensional cycle which can be resolved. The fact that 
a shrinking del Pezzo four-cycle in a CY realizes a theory with exceptional havor symmetry 
leads to a clear candidate for the nature of this new direction. 
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Some remaining issues 


1. Our calculations should extend to the case of real gauge groups, and in particular 
to the theories with extra tensor representations, arising from D3-brane probes of 
F-theory. 

2. “nn” completions of the Seiberg-Witten curve can be seen from the matrix model. 
As explained in [^, including more of the “fractional branes” of the CY singularity 
allows one to determine an embedding of the SW curve in a threefold of a form such 
as 

uv = F{x, y). 

This is important, for example, because it will allow one to identify the resolution 
involved in the extremal transition proposed above. 

3. The held theories obtained on 3-brane probes of F-theory exhibit S-duality. In a 

beautiful series of papers the S-duality of the A/” = 4 theory and its 

M = 1 deformations has been found via the solution of the corresponding matrix 
integral in 0- The symplectic matrix model with the corresponding matter content 
should also have modular behavior in r; it shares the feature with the A/" = 4 theory 
that the path integral over the matter cancels the Vandermonde for the adjoint matrix. 


4. Kazakov [1^] computes “average numbers of holes” and “average lengths of holes” 
in the random surfaces, from the large M solution to the matrix integral. These 
observables exhibit more detailed critical behavior than we have discussed thus far as 
7 , m, Qk are varied. The transition to “torn surfaces” with large holes likely has an 
interpretation in terms of the appearance of a Higgs branch in the gauge theory when 
an rua approaches an eigenvalue of $. It will be interesting to hnd a superpotential 
via which we can £x the moduli at a point on this Higgs branch. 
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